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1.  PROJECT  SUMMARY 

The  objective  of  this  research  program  is  to  study  theoretically  the  underlying  principles  of  solid- 
state  dynamics  and  quantum  mechanical  transport  of  carriers  in  ultrasmall  novel  semiconductor  devices. 
The  areas  of  research  to  be  investigated  are:  1)  theory  of  optical  phonon  modes  in  heterostructures  and 
2)  quantum  transport  in  solids  with  special  emphasis  on  non-perturbative  role  of  high-electric  fields, 
many-body  effects,  and  band-mixing  in  dynamical  processes.  These  problems  will  be  treated  analyti¬ 
cally  through  the  development  of  macroscopic  and  microscopic  physical  models  with  an  emphasis  on 
quantum  mechanical  principles.  Specific  subjects  discussed  in  this  proposal  include  the  effects  of 
confined  and  interface  phonon  modes,  two-electron  quantum  transport  theory  under  hot-electron  condi¬ 
tions,  dielectric  response  function  theory,  and  the  effects  of  band-mixing  in  tunneling.  The  knowledge 
developed  in  this  work  will  be  of  major  importance  in  explaining  the  novel  phenomena  and  fundamen¬ 
tal  questions  relating  to  the  breakdown  of  classical  solid-state  electronics  as  device  dimensional  scales 
are  reduced  to  the  submicron  and  ultra-submicron  regime. 


2 


2.  PROJECT  DESCRIPTION 


2.1  Introduction 

Recent  advances  in  semiconductor  materials  growth  techniques,  such  as  molecular  beam  epitaxy 
(MBE),  metalorganic  chemical  vapor  deposition  (MOCVD),  and  atomic  layer  epitaxy  (ALE),  have 
made  possible  the  fabrication  of  devices  with  one  or  more  dimensions  approaching  the  spacing 
between  planes  of  atoms.  These  nanometer-scale  techniques  have  also  opened  new  possibilities  for 
"band  gap  engineering"  of  novel  semiconductor  devices,  including  heterostructures  with  spatially 
modulated  energy  band  gaps.  By  using  the  variability  of  the  boundary  conditions  which  can  be 
imposed  on  the  wavefunctions,  electrical  and  optical  responses  in  these  structures  can  be  tailored  virtu¬ 
ally  at  will.  As  fabrication  technology  has  allowed  such  ultrasmall  structures  to  be  realized,  many  new 
and  fundamental  questions  have  emerged  concerning  the  underlying  physics  of  small  dimensions  with 
complex,  quantum-scale  boundary  conditions  in  semiconductor  devices.  Important  issues  now  under 
investigation  include  quantum  mechanical  phenomena  such  as  size  quantization,  phase  coherence,  and 
highly  nonequilibrium  transport  in  which  perturbative  treatment  or  the  assumption  of  linear  response  is 
not  applicable  (see,  for  example.  Refs.  1-3).  More  than  ever  before,  it  is  important  that  our  ability  to 
analyze  these  physical  phenomena  in  ultrasmall  devices  keep  in  advance  of  the  technological  develop¬ 
ments. 

In  1990,  the  Office  of  Naval  Research  initiated  sponsorship  of  a  basic  research  program  in  the 
Department  of  Electrical  and  Computer  Engineering  at  North  Carolina  State  University.  The  general 
goal  of  this  research  program  is  to  study  theoretically  solid-state  dynamics  and  quantum  mechanical 
transport  of  carriers  in  ultrasmall  novel  semiconductor  devices.  The  emphasis  is  on  the  development  of 
physical  understanding  of  the  novel  phenomena  through  analytical  approaches  rather  than  numerical 
modeling.  Our  efforts  have  been  devoted  to  investigating  two  important  problems.  The  first  one  is  the 
theory  of  optical  phonon  modes  in  heterostructures.  The  changes  in  longitudinal-optical  (LO)  phonon 
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frequencies,  lifetime,  and  interaction  with  carriers  as  a  result  of  reduced  dimensionality  are  the  main 
subjects  in  this  part  of  research.  Due  to  the  advent  of  lattice  mismatched  strained-layer  (i.e.,  pseu- 
domotphic)  structures,  the  effects  of  strain  on  LO-phonon  modes  have  been  investigated  as  well.  The 
results  will  be  of  major  importance  to  a  wide  variety  of  III-V  semiconductor  devices  in  which  the 
scattering  by  LO-phonon  modes  is  an  important  (at  times,  dominant)  energy  loss  mechanism.  The 
second  topic  is  the  quantum  mechanical  transport  of  charge  carriers  with  specific  emphasis  on  the  role 
of  electron-electron  and  electron-phonon  interactions  in  the  presence  of  high-electric  fields,  many-body 
effects,  and  band-mixing  in  dynamical  processes.  A  novel  formalism  for  treating  Bloch  electron 
dynamics  and  quantum  transport  in  electric  fields  of  arbitrary  strength  and  time  dependence  has  been 
studied  in  an  attempt  to  include  all  quantum  mechanical  effects  collectively  in  the  lowest  order  in  the 
scattering  strength  [1],  Specific  interest  in  this  inelastic  scattering  problem  arises  from  issues  and  ques¬ 
tions  relevant  to  the  role  of  high-electric  fields  in  influencing  transport,  electron  relaxation,  noise  gen¬ 
erations  as  well  as  ionization  processes  in  quantum  wells,  tunnel  barriers,  and  superlattices  [4].  At  the 
same  time,  attempts  are  made  to  study  the  effects  of  quantum  mechanical  principles  on  carrier  transport 
with  more  macroscopic  approaches  such  as  dielectric  function  theory. 

Although  less  than  three  years  old,  this  program  has  already  resulted  in  twenty  refereed  publica¬ 
tions  in  the  literature,  four  additional  manuscripts  are  currently  in  press,  and  five  more  have  been  sub¬ 
mitted  to  major  technical  journals.  Numerous  invited  talks  and  presentations  have  been  and  will  be 
given  at  conferences  and  workshops  throughout  the  United  States  and  in  other  countries.  A  listing  of 
these  publications  is  given  in  APPENDIX  A.  In  addition,  this  research  program  provides  a  perfect 
complement  to  other  aspects  of  our  research  efforts  in  which  numerical  aspects  of  transport  study  arc 
emphasized.  Due  to  the  complicated  geometry  and  boundary  conditions,  a  realistic  answer  for  a  realis¬ 
tic  device  structure  is  attainable  only  through  numerical  approaches.  The  accurate  physical  models  and 
understanding  developed  in  this  research  will  facilitate  the  development  of  the  numerical  models  better 
able  to  explain  experimental  observations  and  to  predict  new  physical  phenomena.  Overall,  this 
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research  program  supported  by  CNR  has  been  efficient  and  productive.  The  quality  of  our  program 
will  continue  to  improve  in  the  future.  In  this  annual  report,  the  progress  and  accomplishments  made 
during  the  past  contract  period  will  be  summarized. 

The  structure  of  the  rest  of  this  report  is  as  follows:  in  Section  2.2,  the  research  results  of  the 
past  contract  period  are  summarized;  Section  2.3  provides  the  list  of  publications  resulting  from  this 
research  supported  by  ONR  during  the  1991  contract  period;  Section  3  contains  background  informa¬ 
tion  describing  project  personnel;  and  the  Appendices  include  a  list  of  refereed  publications  supported 
by  this  research  program  since  1990  and  the  title  page  of  each  paper  published  during  the  1991  con¬ 
tract  period. 

2.2  Summary  of  Research  Results 

This  section  provides  the  current  status  of  our  research  efforts  on  the  theory  of  electron-optical- 
phonon  interactions  and  quantum  mechanical  transport  in  semiconductor  nanostructures.  The  major 
accomplishments  during  the  past  contract  period  are: 

•  We  have  modeled  the  electron-interface-phonon  interaction  with  a  scalar  potential  and  calculated 
the  scattering  strength  in  a  short-period  superlattice  in  an  effort  to  alleviate  controversies  sur¬ 
rounding  the  nature  of  interface  phonon  interaction.  The  results  have  provided  a  criterion  which, 
in  conjunction  with  experiments,  has  been  used  to  examine  and,  subsequently,  to  further  demon¬ 
strate  the  validity  of  the  dielectric  continuum  model; 

•  We  have  studied  the  effects  of  surface  roughness  (i.e.,  non-uniform  thickness)  on  electron- 
optical-phonon  scattering  rates  in  a  rectangular  quantum  wire  and  have  estimated  electron  mobili¬ 
ties  by  applying  an  ensemble  Monte  Carlo  method; 


We  have  investigated  the  electron-interface-phonon  interaction  in  a  double -barrier  heterostructurc 
to  explore  the  possibility  of  achieving  a  simplified  description  for  interaction  Hamiltonians  in 
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multi-heterointerface  structures.  It  is  found  that  when  the  dimensions  of  the  structures  are  larger 
than  approximately  30  A,  the  Hamiltonians  developed  for  a  single  quantum-well  structure  can,  in 
general,  be  applied  to  multi-heterointerface  structures  with  reasonable  accuracy; 

•  We  have  demonstrated  that  optical  phonon  modes  may  be  tailored  through  the  judicious  use  of 
metal-semiconductor  interfaces  in  such  a  way  as  to  dramatically  reduce  unwanted  emission  of 
interface  phonons; 

•  We  have  developed  the  dielectric  response  function  which  properly  incorporates  the  finite  particle 
lifetime  in  heavily-doped  zincblende  semiconductors  at  finite  temperature  and  have  demonstrated 
its  importance  in  electron  transport,  especially,  in  p-type  semiconductors;  and 

•  We  have  derived  the  envelope  functions  for  Bloch  electron  in  the  Wannier  representation  in  the 
presence  of  Slater-Koster  type  impurities  and  have  deduced  the  "transmission  coefficient"  for 
electron  tunneling  through  impurities  and  "band-engineered"  quantum  barriers. 

Each  of  these  results  is  discussed  in  detail  below. 

Role  of  interface-phonon  scattering  in  mesoscale  structures 

Since  the  early  work  of  Fuchs  and  Kliewer  [5],  the  effects  of  confinement  on  LO-phonon  modes 
have  been  studied  theoretically  by  a  number  of  researchers  (including  ourselves).  There  have  been 
indications  that  interactions  by  interface  modes  can  be  significant  and  may  become  the  dominant 
energy  loss  mechanism  in  polar  semiconductors  with  ultra-small  dimensions.  Therefore,  an  appropriate 
treatment  of  electron-interface-phonon  scattering  in  quantized  systems  is  essential  for  the  understanding 
of  election  transport  in  heterostructures.  Recently,  however,  the  role  of  interface  phonons  in  election 
transport  has  been  questioned  [6J.  The  predicted  strong  influence  of  electron-interface-phonon  coupling 
is  obtained  by  applying  the  dielectric  continuum  model  with  electrostatic  boundary  conditions.  There 
is  concern  regarding  the  validity  of  using  the  electrostatic  boundary  conditions  as  well  as  with  describ- 


6 


ing  the  interface  modes  in  terms  of  a  scalar  potential. 

In  order  to  alleviate  controversies  surrounding  this  issue,  we  have  modeled  the  electron-interface- 
phonon  interaction  with  a  scalar  potential  to  calculate  the  ratio  of  electron-interface-phonon  scattering 
to  electron-confined-LO-phonon  scattering  in  a  short-period  superlattice.  The  results  provide  a  criterion 
which,  in  conjunction  with  experiments,  can  be  used  to  examine  the  validity  of  the  dielectric  contin¬ 
uum  model.  The  specific  example  considered  in  this  study  is  a  GaAs/AlAs  short-period  superlattice 
with  the  GaAs  layer  width  of  a  and  the  AlAs  layer  width  of  b.  The  ratio  of  electron-interface-phonon 
scattering  to  electron-confined-LO-phonon  scattering  in  this  superlattice  is  calculated  through  the  Fermi 
golden  rule  by  taking  the  electronic  wavefunction  along  the  growth  direction  (i.e.,  the  z-direction)  as 
that  given  by  Cho  and  Prucnal  for  the  maximum  edge  of  the  first  subband  [7],  To  obtain  an  order-of- 
magnitude  estimate  of  the  ratio  of  electron-interface-phonon  scattering  to  electron-confined-LO-phonon 
scattering,  we  consider  the  ratio  of  principal  factors,  Y,  in  the  integrands  of  the  respective  scattering 
rates.  In  this  calculation,  the  contribution  by  the  S-  mode  is  neglected  due  to  its  relatively  small  con¬ 
tribution  as  compared  with  the  S+  mode.  In  addition,  we  have  included  only  the  lowest  order  confined 
mode  since  it  makes  the  dominant  contribution  among  all  of  the  confined  modes  for  the  electronic 
states  under  consideration.  Figure  1  shows  the  ratio,  Y,  plotted  as  a  function  of  X  (=qa  where  q  is  the 
phonon  momentum  parallel  to  the  heterointerface)  for  a  superlattice  with  b=a  and  b=3a/2.  For  X  <  Jt, 
the  ratio  Y  is  larger  than  unity.  A  principal  factor  contributing  to  the  rapid  increase  in  Y  as  shown  in 
Fig.  1  for  small  values  of  qa  is  the  exponential  increase  in  the  strength  of  the  interface-phonon  mode 
near  the  heterojunction.  In  agreement  with  our  results,  a  surge  in  the  strength  of  the  S+  mode  with 
decreasing  well-width,  a,  is  observed  in  the  recent  low-temperature  Raman  experiment  by  Tsen  et  al. 
[8],  which  appears  to  indicate  the  dominance  of  the  S+  interface-phonon  mode  over  the  confined  LO- 
phonon  modes.  In  addition,  electron  mobilities  measured  by  Zhu  in  short-period  superlattices  show  a 
degradation  as  the  well-width  is  reduced  to  a  very  small  dimension,  which  can  be  attributed  to  an 
increasing  interface-phonon  scattering  as  well  [9].  Therefore,  we  expect  that  electron-interface-phoncn 
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interaction  can  be  described  by  the  dielectric  continuum  model  appropriately  and  should  be  included  in 
the  study  of  low-voltage  (also,  low-temperature  in  general)  operation  characteristic  of  mesoscale  dev¬ 
ices  currently  being  considered  by  the  electron  device  community.  This  point  has  been  proved  further 
by  the  recent  development  in  detailed  microscopic  calculations  of  interface  modes  in  polar  semiconduc¬ 
tors  [10]. 

As  discussed  above,  an  appropriate  treatment  of  the  optical  phonon  modes  in  quantized  systems 
is  essential  for  the  understanding  of  electron  transport  in  heterostructures.  However,  most  of  the 
theoretical  analyses  have  been  confined  to  highly  symmetric  and/or  simple  structures  such  as  single 
quantum  wells  or  superlattices.  Application  of  even  a  simple  macroscopic  model,  not  to  mention  the 
microscopic  ab  initio  models,  becomes  highly  complicated  due  to  the  coupling  between  adjacent  inter¬ 
faces  when  the  structure  has  multiple  heterointerfaces  or  is  asymmetric.  Thus,  it  is  desirable  *o  investi¬ 
gate  the  possibility  of  utilizing  an  approximate  description  of  phonon  modes  and  related  interaction 
Hamiltonians,  while  retaining  the  quantized  nature  of  phonon  modes  in  multi-heterointerface  structures. 
During  the  past  contract  period,  we  have  studied  the  electron-optical-phonon  interaction  in  a 
GaAs/ALAs  double-barrier  heterostructure  to  explore  the  possibility  of  achieving  a  simplified  descrip¬ 
tion  such  as  discussed  above.  The  double-barrier  structure  has  been  chosen  as  a  specific  example  due 
to  its  structural  similarity  (with  the  added  complexity  of  finite  barrier  width)  to  single  quantum  wells 
(i.e.,  infinite  barrier  width)  for  which  quantized  optical  phonon  modes  are  already  well  known.  At  the 
same  time,  this  structure  has  played  a  prominent  role  in  recent  studies  which  demonstrated  the 
existence  and  importance  of  interface-phonon  modes  (see,  for  example,  Ref.  11).  In  our  study,  the 
interaction  Hamiltonians  were  derived  by  applying  the  macroscopic  dielectric  continuum  model  (appli¬ 
cation  of  which  is  justifiable  as  discussed  above)  without  any  further  approximations.  Then,  the  result¬ 
ing  scattering  rates  were  compared  with  those  obtained  by  assuming  the  Hamiltonians  appropriate  for 
single  quantum  wells  for  which  finite  barrier  widths  are  neglected.  Since  the  dielectric  continuum 
model  approximates  the  confined  LO  modes  as  being  completely  confined  to  the  respective  layer  and 
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not  affected  by  or  coupled  to  the  adjacent  layers,  we  have  only  considered  the  influence  of  finite  bar¬ 
rier  widths  on  the  interface  LO  modes  as  well  as  their  interactions  with  electrons.  Figure  2  shows  the 
electron  scattering  rates  via  interface-phonon  emission  as  a  function  of  barrier  thickness,  d2,  when  the 
well  width,  dj,  is  fixed  to  60  A.  Each  of  these  rates  is  evaluated  at  a  value  of  Ey  which  is  chosen  to 
be  just  large  enough  to  permit  a  phonon  emission  so  that  there  is  an  electronic  transition  from  an  initial 
state  to  a  final  state;  i.e.,  for  Fig.  2(a),  Ey  is  fixed  at  37  meV  while,  for  Fig.  2(b),  it  is  37  meV  plus  the 
energy  separation  between  the  first  and  second  subbands.  As  expected,  the  two  rates  merge  smoothly 
with  increasing  d2.  By  the  time  d2  is  approximately  equal  to  d]  (i.e.,  60  A),  the  difference  is  only 
about  20  %.  This  reduces  to  roughly  10  %  with  d2  =  2d]  (=120  A).  In  the  other  limit,  the  difference 
magnifies  rapidly  as  d2  shrinks  below  30  A.  Thus,  the  results  in  Fig.  2  suggest  that  when  d2  >  30  A, 
our  simplified  description  using  Hamiltonians  for  single  quantum  wells  can  be  applied  with  good  accu¬ 
racy  in  place  of  the  more  complicated  expressions  derived  for  each  double-barrier  structure.  When  the 
well  width  is  reduced  from  the  60  A  used  in  Fig.  2,  it  is  expected  that  there  will  be  some  effects  on 
the  coupling  between  the  interface-phonon  modes  in  the  barrier  region.  Our  calculation  shows  that 
when  d]  is  reduced  to  30  A,  only  a  minor  change  is  observed  in  the  phonon  modes  in  the  barrier 
region.  At  the  same  time,  the  electronic  wavefunctions  are  rather  well  confined  in  the  well  region  with 
this  thickness  range.  Hence,  when  the  thicknesses  of  both  the  barrier  and  well  regions  are  larger  than 
30  A  (roughly),  it  is  expected  that  the  electron-interface-phonon  scattering  rates  can  be  well  approxi¬ 
mated  by  the  Hamiltonians  originally  developed  for  single  quantum-well  structures.  This  finding  can 
be  applied  to  multi-heterointerface  structures  in  general  as  well  as  to  the  double-barrier  structures 
investigated  in  this  study.  Application  of  the  simplified  Hamiltonians  to  calculate  the  phonon-assisted 
tunneling  current  in  a  double-barrier  resonant  tunneling  structure  resulted  in  remarkably  good  agree¬ 
ments  with  experimental  data  [11],  which  provides  further  justification  for  the  dielectric  continuum 


model. 


(a) 


(b) 

Figure  2.  Electron  scattering  rates  by  interface-phonon  emission  as  a  function  of  barrier  thickness,  d2 
when  the  well  width,  dh  is  60  A  in  a  GaAs/AlAs  double-barrier  heterostructure  at  300  K.  The  data  in 
(a)  describe  the  intra-subband  scattering  rates  of  the  lowest  resonant  subband  in  the  GaAs  quantum 
well  (i.e.,  1  -»  1),  and  those  in  (b)  show  the  inter-subband  scattering  rates  from  the  lowest  to  the 
second  lowest  subband  (i.e.,  1  -»2).  The  electron  energy  E,(  is  fixed  at  37  meV  for  (a),  while  it  is  37 
meV  plus  the  energy  separation  between  the  first  and  second  subbands  (AE,2)  for  (b).  The  solid  line  is 
obtained  from  the  interaction  Hamiltonians  derived  for  the  double-barrier  structure,  while  the  dashed 
line  is  estimated  by  the  simplified  Hamiltonians  appropriate  for  single  quantum-well  structures. 


The  significance  of  interface  phonon  scattering  with  decreasing  structural  sizes,  as  demonstrated 
above,  poses  a  serious  problem  in  achieving  nanoscale  structures  with  greater  complexity  and  higher 
speed.  Last  year,  we  have  realized  that  elimination  or  dramatic  reduction  of  this  unwanted  energy  loss 
mechanism  may  be  possible  under  a  special  set  of  boundary  conditions.  Based  on  the  dielectric  contin¬ 
uum  model  [5],  it  is  predicted  that  at  the  semiconductor  boundaries  with  an  ideal  metal,  the  optical 
phonon  modes  become  evanescent  and  vanish  unless  they  have  even  polarization  vectors  normal  to  the 
interface.  Thus,  the  phonon  modes  may  be  tailored  through  the  judicious  use  of  metal-semiconductor 
interfaces  in  such  a  way  as  to  drastically  reduce  the  scattering  rates  by  interface  phonons.  This  may 
represent  a  major  breakthrough  for  the  future  generation  electronic  devices.  For  further  demonstration 
of  this  concept,  more  investigation  is  warranted  which  takes  into  account  the  full  dependence  of  pho¬ 
non  suppression  on  the  conductivity  of  the  metal. 

Electron-optical-phonon  interaction  in  realistic  quantum  wires 

The  progress  in  semiconductor  technology  has  provided  the  means  to  fabricate  so-called  quantum 
wires  with  quasi-one-dimensional  (ID)  structures.  It  has  been  suggested  that  quantum  wires  will  exhi¬ 
bit  carrier  mobilities  well  above  106  cm2/V-sec,  but  these  high  values  of  the  mobility  have  not  yet  been 
observed  experimentally.  The  expected  enhancement  of  the  carrier  mobility  in  quantum  wires  should 
stem  from  the  restriction  of  momentum  space  to  one  dimension  as  well  as  the  resulting  reduction  of 
final  states  for  scattered  electrons.  This  point,  however,  needs  to  be  clarified  since  electron  interactions 
with  the  dominant  optical  phonon  modes  are  strongly  modified  by  phonon  confinement  and  localization 
as  well. 

During  the  1990  contract  period,  we  have  reported  the  derivation  of  electron-optical-phonon 
interaction  Hamiltonians  in  a  rectangular  quantum  wire,  in  which  the  phonon  confinement  and  localiza¬ 
tion  (i.e.,  confined  and  interface  modes)  are  considered  by  applying  the  dielectric  continuum  model 
112).  In  this  study,  an  "ideal"  structure  is  assume  which  means  that  the  surface  roughness  of  the 
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quantum  wire  is  neglected  (i.e.,  a  constant  thickness  along  the  structure).  Real  quantum  wires,  how¬ 
ever,  always  have  variable  thickness  along  the  structure.  Generally,  electron  beam  or  X-ray  lithogra¬ 
phy  on  a  quantum  well  surface,  with  subsequent  reactive  ion  (beam)  etching,  is  used  for  quantum  wire 
fabrication.  State-of-the-art  technology  allows  the  quantum  well  width  to  be  controlled  down  to  one 
monolayer.  However,  the  etching  results  in  quantum  wires  with  thicknesses  varying  within  several 
percent  along  the  structure.  Even  the  most  sophisticated  technologies  fail  to  provide  abrupt  interfaces 
between  the  layers.  Thus,  it  appears  that  all  current  and  foreseeable  future  technologies  for  fabricating 
quantum  wires  do  not  assure  the  possibility  to  creating  ideal  structures  with  constant  thickness.  The 
variation  of  quantum  wire  thickness  results  in  the  variation  of  subband  energies  and,  consequently, 
election  scattering  is  no  longer  energetically  coherent  in  different  parts  of  a  quantum  wire.  This  may 
significantly  affect  the  unique  features  of  ideal  ID  systems  (such  as  diverging  density  of  states  at  the 
subband  minimum)  and  alter  the  carrier  transport  properties  in  realistic  situations. 

As  a  generalization  of  our  previous  results,  we  have  started  an  investigation  on  the  accurate 
description  of  election  transport  in  real  quantum  wire  structures.  For  simplicity,  the  real  quantum 
wires  are  assumed  to  have  a  smooth,  harmonic  variation  of  thickness  along  the  structure  so  that  the 
characteristic  length  of  the  fluctuations  is  much  greater  than  the  de  Broglie  wavelength  AB.  The  varia¬ 
tion  amplitudes  5Ly  and  5LZ  have  been  chosen  as  a  fraction  ranging  from  1  %  to  10  %  of  the 
corresponding  thicknesses  Ly  and  Lz.  We  have  considered  three  cases  with  different  mutual  phases  of 
the  variation  of  Ly  and  Lz:  i.e.,  Lz  =  Lm  +  8Lzcos(Kx),  Lz  =  Lzo  +  5Lzsin(Kx),  and  Lz  = 
L m  +  8Lzcos(2Kx)  with  Ly  =  Ly0  +  5LyCOs(Kx)  in  all  three  cases.  The  scattering  rates  k  by  confined 
and  interface  optical  phonon  modes  are  calculated  with  a  give  Ly  and  Lz  based  on  the  ID  interaction 
Hamiltonians  developed  previously  [12],  and  properly  averaged  over  the  quantum  wire  length  L  with 
varying  thicknesses: 


11 


L 

XO^L^L^)  =  1  J  dx  X[Ex(x),Ly(x),Lz(x)]  ,  (1) 

L  o 

where  Ej  =  E0  +  5E(x),  E0  =  E  -  <Eji>  is  the  electron  kinetic  energy  corresponding  to  a  given  total 
energy  E  and  average  subband  position  <Ejj>,  5E(x)  is  the  electron  kinetic  energy  variation  (opposite  to 
the  variation  of  subband  position)  due  to  the  variation  of  thickness.  It  has  been  found  that  the  results 
do  not  show  significant  differences  between  these  models  as  long  as  it  is  in  the  limit  of  smooth  fluctua¬ 
tions  (i.e.,  KAb  c  1).  A  typical  electron  transition  rate  from  the  first  subband  to  elsewhere  due  to 
emission  of  confined  LO  phonons  is  presented  in  Fig.  3  for  ideal  and  real  quantum  wires.  A  GaAs 
quantum  wire  embedded  in  AlAs  is  assumed  with  the  dimensions  of  Ly=150  A  and  Lz=250  A.  The 
forward  and  backward  scatterings  are  plotted  as  separate  scattering  mechanisms  in  order  to  reveal  the 
polar  character  of  electron-phonon  interaction.  The  scattering  rates  by  interface  modes  are  not  con¬ 
sidered  since  they  are  more  than  an  order  of  magnitude  smaller  than  those  by  confined  LO  modes,  as 
discussed  in  our  previous  report.  For  a  reasonable  5  %  variation  in  amplitude  (i.e.,  8Ly/Ly0  = 
SLJLjo  =  0.05),  one  can  see  in  Fig.  3  that  the  scattering  rates  are  considerably  modified  in  real  quan¬ 
tum  wires  compared  to  the  ideal  wire.  It  is  interesting  to  note  that  the  well-pronounced  resonant  nature 
of  electron  scattering  in  the  form  of  multiple  sharp  peaks  (which  is  one  of  the  unique  features  in  an 
ideal  ID  system)  is  significantly  broadened  and  reduced  (for  the  first  peak)  or  even  completely  disap¬ 
peared  (for  the  higher  peaks).  Similar  curves  with  a  different  degree  of  broadening  are  obtained  for 
other  variation  amplitudes.  Even  a  variation  in  amplitude  as  small  as  1  %  leads  to  the  disappearance 
of  the  divergence  at  the  resonant  energies  and  a  dramatic  smearing  out  of  the  peak -like  structure  on  the 
curves  discussed. 

The  calculated  scattering  rates  for  an  ideal  quantum  wire  have  been  included  in  the  Monte  Carlo 
program  and  electron  transport  has  been  simulated  at  various  electric  fields  as  well.  The  aim  of  these 
simulations  is  to  reveal  the  role  of  electron  intersubband  scattering  due  to  confined  and  interface  pho- 
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Figure  3.  Electron  scattering  rate  due  to  LO  phonon  emission  in/ffom  the  first  subband  as  a  function 
of  the  electron  kinetic  energy.  Dashed  curves  are  for  the  ideal  quantum  wire;  and  solid  curves  for  the 
real  quantum  wire  with  5  %  thickness  variation  along  the  structure.  Upper  curves  represent  forward 
electron  scattering  (wavevector  direction  does  not  change);  and  lower  curves  represent  backward 
scattering  (wavevector  change  its  direction  to  opposite).  A  GaAs  quantum  wire  embedded  in  AlAs  is 
assumed  with  the  dimensions  of  Ly=150  A  and  1^=250  A  at  300  K. 


Figure  4.  Electron  drift  velocity  as  a  function  of  electric  field.  The  data  points  represent  the  results 
obtained  for  the  quantum  wire  and  the  curves  are  those  for  bulk  GaAs.  Upper  dots  and  curves  are  at 
77  K  and  the  lower  dots  and  curves  are  at  300  K. 
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nons  in  electron  transport  phenomena  in  quantum  wires,  and  to  examine  the  feasibility  of  the  suggested 
high  mobilities  in  quantum  wire  structures.  Under  the  conditions  of  intensive  intersubband  scattering 
at  300  K,  electrons  do  not  exhibit  the  expected  high  mobilities,  as  one  can  see  from  Fig.  4  where 
velocity-field  dependences  for  the  quantum  wire  and  bulk  undoped  GaAs  are  plotted  for  comparison. 
At  this  temperature,  electron  scattering  by  confined  LO  phonons  determines  electron  low-field  mobility. 
The  role  of  acoustic  and  impurity  scattering  as  well  as  scattering  by  interface  phonons  is  negligible. 
Electron  mobility  enhancement  in  ideal  quantum  wires  is  seen  in  Fig.  4  at  77  K  and  low  electric  fields; 
in  that  case  electrons  do  not  undergo  effective  scattering  in  the  "passive”  region  below  the  optical  pho¬ 
non  energy.  However,  at  77  K  the  probability  of  absorption  of  optical  phonons  is  of  the  same  order  as 
that  of  acoustic  scattering  and  the  electron  mobility  is  controlled  not  only  by  the  electron-optical- 
phonon  scattering  but  also  by  the  acoustic  phonon  scattering  which  is  not  considered  in  the  present 
work.  Therefore,  we  observe  higher  mobility  than  it  should  be.  It  must  be  noted  that  electron  scatter¬ 
ing  by  interface  phonons  does  not  influence  considerably  electron  steady-state  high-field  transport 
However,  interface  phonons  which  have  energies  different  from  those  of  confined  LO  phonons  yield 
very  fast  energy  dissipation  and  are  essential  to  study  the  relaxation  of  the  electron  energy  distribution 
function  at  low  electric  fields.  In  the  case  when  interface  phonon  scattering  is  not  considered,  the 
injected  electron  distribution  function  exhibits  sharp  peaks  at  energies  which  are  multiples  of  the 
confined  LO  phonon  energy  and  does  not  relax  to  the  equilibrium  distribution  function  in  the  absence 
of  an  electric  field. 

Effects  of  finite  particle  lifetime  on  electron  transport  in  heavily-doped  semiconductors 

During  the  past  contract  period,  we  have  initiated  a  study  on  the  dielectric  response  function 
theory  to  advance  our  understanding  of  many-body  effects,  such  as  plasmon-phonon  coupling  and 
dynamic  screening,  within  the  frame  work  of  the  random  phase  approximation  (RPA)  or  the  self- 
consistent  field  approach  (SCF).  It  is  well  known  that  the  dynamic  response  of  a  free  electron  gas  is 
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described  by  the  Lindhard  dielectric  function  in  the  RPA  or  the  SCF  when  the  electron  lifetime  x  is 
infinite.  The  free  electron  gas  formula  is  easily  applicable  to  heavily-doped  n-type  semiconductors 
since  they  have  a  single  conduction  band  with  s-like  symmetry.  For  p-type  semiconductors,  the  situa¬ 
tion  is  much  more  complicated  because  of  the  p3/2-like  symmetry  of  the  valence  bands.  Recently, 
Yevic  and  Bardyszewski  [13,14]  evaluated  the  Lindhard  dielectric  function  for  a  p-type  semiconductor 
at  finite  temperature.  These  authors  showed  that  the  structure  of  the  spectral  density  function  is  quite 
different  from  that  for  zero  temperature.  However,  they  still  neglected  the  effects  of  the  finite  hole 
lifetime  even  for  finite  temperature.  Although  the  Lindhard  dielectric  function  inherently  incorporates 
the  effects  of  single  particle  excitation  (Landau  damping)  including  both  intra-  and  interband  transi¬ 
tions,  it  does  not  take  into  account  the  finiteness  of  particle  lifetime  due  to  the  collision  effects  such  as 
impurity  scattering  and  phonon  scattering.  At  finite  temperature,  these  effects  are  very  important  since 
such  scattering  disturbs  the  collective  motion  of  quasi -particles  and  therefore  the  dynamic  response 
function  will  be  affected  significantly.  As  a  result,  electron  transport  in  such  environments  will  be 
influenced  strongly  as  well.  A  detailed  understanding  of  the  behavior  of  electrons  in  heavily-doped 
semiconductors  is  increasingly  important,  especially  for  the  design  of  high  speed  HBTs  and  HETs, 
since  electron  transport  in  the  base  region  is  one  of  the  critical  factors  limiting  the  maximum  switching 
speed  of  such  devices. 


By  modifying  an  approach  proposed  by  Mermin  for  free  electron  gas  [15],  we  have  developed 
the  dielectric  response  function  which  property  incorporates  the  finite  particle  lifetime  in  heavily-doped 
zincblende  semiconductors  (denoted  as  the  Lindhard-Mermin  dielectric  function  hereinafter).  The 
dielectric  response  function  for  a  polar  semiconductor  comprises  three  terms  which  can  be  written  as 


e(q,co)  =  £«,  +  ellU(q,co)  -  p(q,w)  . 


The  first  term  is  the  high  frequency  dielectric  constant  with  vanishingly  small  free  carrier  density, 
which  arises  from  the  polarization  of  valence  electrons  due  to  virtual  interband  transitions.  The  second 


14 


term  is  the  lattice  dielectric  function  usually  described  by  the  well-known  Lyddane-Sachs-Teller  rela¬ 
tionship: 


eUtt(q,w)  =  e» 


<°lo  ~ 

2  2~ 
o>ro  ~  o> 


(3) 


Here,  the  spatial  dispersion  and  damping  of  phonons  due  to  multiple  phonon  processes  are  neglected  to 
avoid  unnecessary  details.  The  third  term  in  Eq.  (2)  is  attributed  to  the  contribution  of  the  carrier 
polarization  where  p(q,co)  represents  the  polarization.  When  the  lifetime  of  the  carrier  gas  is  infinite, 
the  first  and  third  terms  give  the  well-known  Lindhard  dielectric  function: 


eL(q,co)  =  e^~ 


q  lUk 


fn(k)  -  f](k+q) 

E„(k)  -  E,(k  +  q)  -  hco  j 


(4) 


where  E„  and  fn  are  the  particle  energy  and  the  Fermi  distribution  in  the  n-th  band.  The  summation  on 
n  and  1  indices  correspond  to  just  n=l=c  (i.e.,  the  conduction  band)  for  n-type  semiconductors.  For  p- 
type  semiconductors,  the  indices  correspond  to  either  the  heavy  hole  (hh)  or  the  light  hole  (lh)  band. 
The  overlap  factor  IF*jk+ql2  which  is  obtained  from  the  lowest  kp  terms  is  1  for  n  =  1  =  c;  1  - 
3q2sin20/4 1  k+q  1 2  for  n  =  1  =  hh.Ih;  and  Bq^in2©#  I  k+q  1 2  for  n,l  =  hh.lh  with  n*l.  When  the  finite 
lifetime  x  is  considered,  the  third  term  in  Eq.  (2)  is  modified  to  give 

lm.  ,  (l+r7on)[eL(q,co+iVT)  -  ej 

e^Cq.co)  =  e*.  + - - - = -  .  (5) 

1  +  07oK)[eL(q.G>+i/T)  -  eJ/[eL(q,0)  -  ej 


The  evaluation  of  the  appropriate  value  of  the  relaxation  time  x  is  quite  complicated  when  considering 
the  validity  of  relaxation  time  approximation.  However,  since  x  is  related  to  the  conductivity,  an 
approximate  value  can  be  estimated  from  experimental  data  for  mobility.  Once  Eq.  (5)  is  obtained,  the 
total  inelastic  scattering  (by  coupled  plasmon-phonon  interactions)  rate  is  given  by 
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for  an  electron  with  energy  E  in  the  conduction  band,  where  nu  is  the  Bose  function  for  the  phonon 
distribution. 

As  a  specific  example,  we  have  chosen  p-type  GaAs  to  demonstrate  the  significance  of  finite  x  in 
the  electron  transport  studies.  The  spectral  density  functions  are  obtained  for  three  different  values  of 
ft /x  (0,  0.02  eV,  and  0.05  eV)  and  a  hole  density  of  p  =  2xl018  cm-3  at  300  K.  For  these  conditions, 
the  hole  mobility  is  about  100  cm-2/V-sec,  which  corresponds  to  ft/x=0.023  eV.  For  small  values  of  q 
(-1x10s  cm-1),  the  plasmon  damping  is  small  when  ft/x=0,  as  shown  previously  [13].  This  is  due  to 
the  smallness  of  the  Landau  damping.  However,  as  expected,  when  the  relaxation  time  is  set  to  a  finite 
value,  the  plasmon  mode  is  heavily  damped  and  broadened.  For  ft/x=0.05  eV,  the  plasmon  feature 
almost  completely  disappears.  We  also  noticed  that  the  anti-screening  effect  is  weakened  and  the 
phonon-like  mode  is  also  damped  at  a  small  q.  As  q  becomes  larger  (>lxl06  cm-1),  collision  effects 
(due  to  the  finite  hole  lifetime)  become  less  dominant  and  Landau  damping  dominates.  Thus,  for 
minority  electron  transport,  the  finite  hole  lifetime  influences  the  long  range  interaction  (i.e.,  small  q) 
between  an  electron  and  a  hole  gas  significantly,  while  the  short  range  interaction  is  not  affected  appre¬ 
ciably.  The  effects  of  finite  hole  lifetime  on  the  inelastic  scattering  rates  of  minority  electrons  calcu¬ 
lated  from  Eq.  (6)  are  shown  in  Fig.  5.  The  collision  effects  in  the  valence  bands  clearly  enhance  the 
transition  rates  of  minority  electrons.  This  result  is  attributed  mainly  to  the  larger  value  of  the  spectral 
density  function  at  small  q  and  co  due  to  the  collision-induced  broadening,  as  discussed  above. 
Although,  only  the  case  of  p  =  2xl018  cm-3  at  300  K  is  reported  here,  the  collision  effects  are  impor¬ 
tant  in  most  cases  for  heavily-doped  semiconductors  (in  the  n-type  cases,  the  influence  is  generally 
smaller  but  still  sizable).  As  the  doping  density  increases,  the  Landau  damping  and  collision  effects 
(i.e.,  finite  particle  lifetime)  both  become  larger  at  finite  temperature.  In  this  sense,  the  simple  plasmon 
model,  which  has  been  used  in  some  Monte  Carlo  simulations  is  highly  questionable. 
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Figure  5.  Inelastic  scattering  rates  of  minority  electrons  as  a  function  of  electron  energy  for 
p  m  2xl018  cm-3  at  300  K.  Note  that  the  scattering  rates  are  not  zero  even  for  E  =  0  due  to  the  ther¬ 
mal  broadening. 
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Electric  field  assisted  quantum  transport  in  solids 

When  we  deal  with  carrier  transport  in  semiconductor  devices  with  ultrasmall  dimensions,  it  is 
necessary  to  study  quantum  effects  of  carriers  in  semiconductors  subjected  to  rapidly  varying,  spatially 
inhomogeneous  electric  fields  and  non-steady  state  temporal  conditions.  At  this  scale,  many  new  and 
interesting  effects,  which  are  currently  treated  as  only  high  order  corrections,  will  begin  to  dominate 
the  device  performance  and  the  semi -classical  Boltzmann  equation  approach  becomes  highly  question¬ 
able. 

A  novel  formalism  for  treating  Bloch  electron  dynamics  and  quantum  transport  in  inhomogeneous 
electric  fields  of  arbitrary  strength  and  time  dependence  was  developed  previously  by  Iafrate  et  al. 
[1,16].  In  this  formalism,  the  electric  field  is  described  through  the  use  of  vector  potential.  This 
choice  of  gauge  leads  to  a  natural  set  of  basis  functions  for  describing  Bloch  electron  dynamics;  in 
addition,  a  basis  set  of  localized,  electric  field  dependent  Wannier  functions  is  established  and  utilized 
to  derive  a  quantum  "Boltzmann  equation"  which  includes  explicit  band-mixing  transients  such  as 
effective  mass  dressing  and  Zener  tunneling.  In  this  localized  Wannier  representation,  the  differential 
equation  for  the  envelope  function  fn(r,t)  in  the  Wannier  representation  is  given  by 

^  a^t)  =  en(-/V-_LAo)  fn(r ,t)  -  eE0-2;  RnnfiV-lA0)  fn-(r,t)  (7) 

+  IZeVm<r,r)fn<r,t), 

n'  1' 

where  e„(K)  is  the  n-th  Bloch  energy  band  with  crystal  momentum  K,  Rnn'  are  the  band-mixing  terms 
terms  due  to  the  homogeneous  electric  field,  are  the  matrix  elements  of  the  inhomogeneous  poten¬ 
tial  with  respect  to  the  localized  basis,  and  A0  is  the  vector  potential  due  to  the  field  of  spatially  homo¬ 
geneous  part.  In  this  representation,  the  complete  set  of  Wannier  functions  is  used  as  basis  functions. 
These  functions  have  the  desirable  property  of  being  localized  which  makes  them  most  convenient  for 
discussing  a  localized  perturbation  and  the  localized  states  arising  from  this  perturbation. 
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As  an  application  of  the  formalism,  we  have  studied  the  solution  to  Eq.  (7)  for  the  case  of  a 
nearest  neighbor  tight-binding  band  at  zero  homogeneous  field  with  Slater-Koster  impurities  as  inho¬ 
mogeneous  fields.  In  this  case,  the  time-independent  envelope  function  gn(r)  can  be  written  as 

g„(r)  =  e'ur  +  ££  Vn(l',l)  Wn(r,l)  gn(l') , 

i  r  w 


where  Wn(r,l)  is  the  Green’s  function.  The  Green’s  function  has  been  derived  within  the  nearest 
neighbor  tight-binding  approximation.  The  impurity  energy  levels  within  the  energy  bandgap  as  well 
as  the  "transmission  coefficients"  for  a  Bloch  electron  tunneling  through  one  and  two  Slater-Koster 
type  impurities  were  obtained  respectively.  It  is  found  that  the  transmission  coefficient  T  is  always 
less  than  one  in  the  one  impurity  case.  However,  resonant  tunneling  occurs  in  the  two  Slater-Koster 
impurity  case.  The  transmission  coefficient  for  an  electron  to  tunnel  through  two  impurities  with 
impurity  strength  Vc  and  V,  is  found  to  be 


T  = 


_ I _ 

1  +  [2absin(ul)  -  (a+b)cos(ul)]2  +  (a-b)2sin2(ul) 


(9) 


where  a  and  b  are  band  structure  dependent  parameters,  u  is  the  momentum  of  the  electron,  and  1  is  the 
distance  between  two  impurities.  For  the  nearest  neighbor  tight-binding  band,  ^(k^^EjCOsCka); 
a=V0/[2e1sin(u)];  and  b=V)/[2eisin(u)].  For  the  impurities  of  identical  strength, 


T  = 


1  +4[ 


-]2[: 


2e1sin(u)  2eisin(u) 


sin(ul)  -  cos(ul)]2 


(10) 


The  resonant  transmission  (T=l)  occurs  when 


Vq 

2E!Sin(u) 


sin(ul)  -  cos(ul)  =  0  . 


(11) 


Clearly  at  energies  determined  by  Eq.  (11),  electrons  resonantly  tunnel  through  impurity-based  double 


barriers. 


By  treating  a  different  crystal  (lattice-matched)  as  many  impurities,  we  have  examined  the  tran¬ 
sport  of  electron  through  a  double  heterojunction  structure  (i.e.,  band-engineered  quantum  barrier  or 
well)  and  a  single  heterojunction  (i.e.,  infinite  impurities)  as  well.  The  transmission  coefficient  T  for 
an  electron  to  tunnel  through  these  heterostructures  have  been  derived.  For  a  double  heterojunction 
structure  (described  as  L+l  impurities), 

where  c  and  d  are  determined  by  the  momentum  of  electron  u  and  u'  in  the  two  different  crystals, 
respectively.  With  some  assumptions  for  the  energy  bands,  the  transmission  coefficient  T  is  simplified 
as 


T  = 


1 


1  + 


sinV(L4-l)] 

sii.2(u')sin2(u) 


(cos(u')  -  cos(u)]2 


_ 1 _ 

j  +  ,  Vq  )2  sin2[uU+l)] 
2ei  sin2(u')sin2(u) 


03) 


It  predicts  the  resonant  tunneling  when 


sin[u'(L+l))  =  0  u'  =  -j— ,  (14) 

where  m  is  an  integer. 

This  formalism  has  also  been  applied  to  quantum  transport  and  Bloch  oscillations  in  a  homogene¬ 
ous  electric  field.  In  the  case  where  no  inhomogeneity  is  present,  it  has  been  found  that  in  a  one-band 
model  the  homogeneous  electric  field  leads  to  inherent  localization  of  charge  in  periodic  potential,  and 
the  energy  forms  a  Staric  ladder  spectrum.  This  field-induced  localization  is  explicitly  revealed  by  the 
time  evolution  of  the  envelope  function 

f„0,t)  =  I  K(l-r,t-t0  fn(i',tO  =  £  e*®-™  J,-K$)  fnG'.0) ,  (15) 


where  the  propagator  KO-r.t-t'Jse'^-1  4)Jh<^),  with  ^=(W/hcoB)sin(coBt/2)  and  the  Bessel  function 
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is  a  localized  function  for  small  £  (or  high  electric  field).  Here,  W  is  the  band  width  and  coB  is 
the  Bloch  frequency.  Equation  (15)  predicts  localization  of  electron  within  a  single  lattice  site  when 
W/1kob  <1.  In  a  three  dimensional  solid,  in  which  the  lattice  spacing  a  is  on  the  order  of  an  A  and 
the  band  width  W  is  on  the  order  of  an  eV,  the  confinement  of  an  electron  in  a  lattice  site  requires  a 
very  high  electric  field  strength  (-107  V/cm).  On  the  other  hand,  in  a  supcrlattice  with  a  period  of 
-100  A  and  W  in  the  range  of  0.01-0.1  eV,  the  electronic  motion  could  be  restricted  to  a  distance 
comparable  to  the  superlattice  period  for  moderately  high  electric  fields  (-105  V/cm).  Localization  in  a 
superlattice  under  such  a  field  has  been  observed  by  Mendez  et  al.  [17J.  At  such  a  high  electric  field, 
the  possible  path  for  electron  to  escape  from  the  lattice  site  will  be  the  transition  to  another  band  (or 
Zener  tunneling). 

2.3  Publications  and  Presentations 

During  the  last  year,  this  program  has  resulted  in  eleven  refereed  publications  in  the  literature, 
four  additional  manuscripts  are  currently  in  press,  and  five  more  have  been  submitted  to  major  techni¬ 
cal  journals.  Eight  presentations  and  invited  talks  have  been  given  at  conferences  and  workshops.  The 
following  paragraphs  summarize  the  publications  and  presentations  made  under  this  program  during  the 
past  contract  period. 

A.  Refereed  Publications 

M.  A.  Stroscio,  K.  W.  Kim,  G.  J.  Iafrate,  M.  Dutta,  and  H.  L.  Grubin,  "Reduction  and  Control  of  Ine¬ 
lastic  Longitudinal-Optical  Phonon  Scattering  in  Nanoscale  and  Mesoscopic  Device  Structures," 
Proceedings  of  the  1991  lnt.  Semiconductor  Device  Research  Symp.,  pp.  87-91,  1991. 

T.  Kaneto,  K.  W.  Kim,  and  M.  A.  Littlejohn,  "Effects  of  Finite  Hole  Lifetime  on  Electron  Transport  in 
P-Type  Zincblende  Semiconductors,"  Proceedings  of  the  1991  lnt.  Semiconductor  Device  Research 
Symp.,  pp.  235-238,  1991. 

V.  V.  Mitin,  R.  Mickevicius,  M.  A.  Stroscio,  G.  J.  Iafrate,  and  K.  W.  Kim,  "Electron  Dynamics  in 
Quantum  Wires,”  Proceedings  of  the  1991  lnt.  Semiconductor  Device  Research  Symp.,  pp.  111-114, 
1991. 

M.  A.  Stroscio,  G.  J.  Iafrate,  K.  W.  Kim,  M.  A.  Littlejohn,  H.  L.  Grubin,  V.  V.  Mitin,  and  R. 
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Mickevicius,  "Role  of  Phonon  Confinement  in  Nanoscale  Systems,"  in  Nanostructures  and  Mesoscopic 
Systems,  edited  by  W.  P.  Kirk  and  M.  A.  Reed  (Academic  Press,  San  Diego,  Calif.,  1992),  Chap.  8, 
pp.  379-386. 

K.  W.  Kim,  M.  A.  Littlejohn,  M.  A.  Stroscio,  and  G.  J.  Iafrate,  "Transition  from  LO-Phonon  to  SO- 
Phonon  Scattering  in  Mesoscale  Structures,"  Semicond.  Sci.  Technol.  7,  B49  (1992). 

R.  Mickevicius,  V.  V.  Mitin,  K.  W.  Kim,  and  M.  A.  Stroscio,  "Electron  High-Field  Transport  in 
Multi-Subband  Quantum  Wire  Structures,"  Semicond.  Sci.  Technol.  7,  B299  (1992). 

S.  Das  Saima,  V.  B.  Campos,  M.  A.  Stroscio,  and  K.  W.  Kim,  "Confined  Phonon  Modes  and  Hot 
Electron  Energy  Relaxation  in  Semiconductor  Microstructures,"  Semicond.  Sci.  Technol.  7,  B60 
(1992). 

M.  A.  Stroscio,  K.  W.  Kim,  G.  J.  Iafrate,  M.  Dutta,  and  H.  L.  Grubin,  "Dramatic  Reduction  of  the 
Longitudinal-Optical  Phonon  Emission  Rate  in  Polar- Semiconductor  Quantum  Wires,"  Phil.  Mag.  Lett. 
65,  173  (1992). 

R.  Mickevicius,  V.  V.  Mitin,  K.  W.  Kim,  and  M.  A.  Stroscio,  "Electron  Intersubband  Scattering  in 
Real  Quantum  Wires,"  Superlatt.  Microstruct.  11,  277  (1992). 

R.  Mickevicius,  V.  V.  Mitin,  K.  W.  Kim,  M.  A.  Stroscio,  and  G.  J.  Iafrate,  "Electron  Intersubband 
Scattering  by  Confined  and  Localized  Phonons  in  Real  Quantum  Wires,"  J.  Phys.:  Condens.  Matter  4, 
4959  (1992). 

K.  W.  Kim,  A.  R.  Bhatt,  M.  A.  Stroscio,  P.  J.  Turley,  and  S.  W.  Teitsworth,  "Effects  of  Interface  Pho¬ 
non  Scattering  in  Multi-Heterointerface  Structures,"  J.  Appl.  Phys.  72,  2282  (1992). 

H.  Qiang,  F.  H.  Poliak,  C.  M.  Sotomayor-Torres,  W.  Leitch,  A.  H.  Kean,  M.  A.  Stroscio,  G.  J.  Iafrate, 
and  K.  W.  Kim,  "Size  Dependence  of  the  Electron-Optical  Phonon  Coupling  in  GaAs/Ga07Alo.3As 
Single  Quantum  Wells,"  accepted  for  publication  in  Appl.  Phys.  Lett. 

T.  Kaneto,  K.  W.  Kim,  and  M.  A.  Littlejohn,  "Dielectric  Response  Functions  of  heavily-Doped  Zinc- 
blende  Semiconductors  with  Finite  Particle  Lifetime,"  accepted  for  publication  in  J.  Appl.  Phys. 

M.  A.  Stroscio,  K.  W.  Kim,  G.  J.  Iafrate,  M.  Dutta,  and  H.  L.  Grubin,  "Reduction  of  Inelastic 
Longitudinal-Optical  Phonon  Scattering  in  Narrow  Polar-Semiconductor  Quantum  Wells,"  accepted  for 
publication  in  Proc.  SPIE  1675  (1992). 

T.  Kaneto,  K.  W.  Kim,  and  M.  A.  Littlejohn,  "Minority  Carrier  Transport  in  Heavily-Doped  p-Type 
Semiconductors,"  submitted  to  Phys.  Rev.  B. 

A.  R.  Bhatt,  K.  W.  Kim,  M.  A.  Stroscio,  G.  J.  Iafrate,  M.  Dutta,  H.  L.  Grubin,  R.  Haque,  and  X.  T. 
Zhu,  "Reduction  of  Interface  Phonon  Modes  Using  Metal-Semiconductor  Heterostructures,"  submitted 
to  J.  Appl.  Phys. 

M.  U.  Erdogan,  K.  W.  Kim,  and  M.  A.  Stroscio,  "Effects  of  Band  Mixing  on  Hole  Tunneling  Times  in 
GaAs/AlAs  Double-Barrier  Heterostructures,"  submitted  to  Appl.  Phys.  Lett. 
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T.  Kaneto,  K.  W.  Kim,  and  M.  A.  Littlejohn,  "A  Comparison  of  Minority  Electron  Transport  in 
Ino.s3Gao.47As  and  GaAs,"  submitted  to  Appl.  Phys.  Lett. 

M.  A.  Stroscio,  G.  J.  Iafrate,  K.  W.  Kim,  A.  R.  Bhatt,  M.  Dutta,  and  H.  L.  Grubin,  "Reduction  in 
Longitudinal-Optical  Phonon  Emission  Rate  in  Polar-Semiconductor  Quantum  Wires  and  Quantum 
Wells,”  to  be  published  in  the  Proceedings  of  the  7th  Inti.  Conf.  on  Phonon  Scattering  in  Condensed 
Matter  (August,  1992,  Ithaca,  New  York). 

J.  He  and  G.  J.  Iafrate,  "Multi-Band  Theory  of  Zener  Tunneling  in  Solids,"  submitted  to  Phys.  Rev.  B. 


B.  Conference  Presentations 

M.  A.  Stroscio,  K.  W.  Kim,  G.  J.  Iafrate,  M.  Dutta,  and  H.  L.  Grubin,  "Reduction  and  Control  of  Ine¬ 
lastic  Longitudinal-Optical  Phonon  Scattering  in  Nanoscale  and  Mesoscopic  Device  Structures 
(Invited),"  presented  at  the  1991  Inti.  Semiconductor  Device  Research  Symp.  (December,  1991,  Char¬ 
lottesville,  Virginia). 

T.  Kaneto,  K.  W.  Kim,  and  M.  A.  Littlejohn,  "Effects  of  Finite  Hole  Lifetime  on  Electron  Transport  in 
P-Type  Zincblende  Semiconductors,"  presented  at  the  1991  Inti.  Scmiconuuctor  Device  Research 
Symp.  (December,  1991,  Charlottesville,  Virginia). 

V.  V.  Mitin,  R.  Mickevicius,  M.  A.  Stroscio,  G.  J.  Iafrate,  and  K.  W.  Kim,  "Electron  Dynamics  in 
Quantum  Wires  (Invited),"  presented  at  the  1991  Inti.  Semiconductor  Device  Research  Symp. 
(December,  1991,  Charlottesville,  Virginia). 

K.  W.  Kim,  A.  Bhatt,  H.  Grubin,  T.  Zhu,  M.  Dutta,  R.  Haque,  G.  J.  Iafrate,  and  M.  A.  Stroscio, 
"Reduction  of  Interface  Phonon  Scattering  in  Polar  Mesoscopic  Systems,"  presented  at  the  March 
Meeting  of  the  American  Physical  Society  (March,  1992,  Indianapolis,  Indiana),  Bull.  Am.  Phys.  Soc. 
37,  397  (1992). 

M.  A.  Stroscio,  K.  W.  Kim,  G.  J.  Iafrate,  M.  Dutta,  and  H.  L.  Grubin,  "Reduction  of  Inelastic 
Longitudinal-Optical  Phonon  Scattering  in  Narrow  Polar-Semiconductor  Quantum  Wells,”  presented  at 
the  SPIE  Symp.  on  Compound  Semiconductor  Physics  and  Devices  (March,  1992,  Somerset,  Ne  v  Jer¬ 
sey). 

M.  A.  Stroscio,  G.  J.  Iafrate,  K.  W.  Kim,  A.  Bhatt,  M.  Dutta,  and  H.  L.  Grubin,  "Reduction  of 
Carrier-Interface-Phonon  Scattering  Near  Metal-Semiconductor  Interfaces  for  Phonon  Frequencies 
Below  the  Plasmon  Frequencies,"  presented  at  the  19th  IEEE  Inti.  Conf.  on  Plasma  Science  (June, 
1992,  Tampa,  Florida). 

X.  T.  Zhu,  H.  Goronkin,  G.  Maracas,  R.  Droopad,  M.  Stroscio,  G.  Ialrate,  K.  Kim,  and  A.  Bhatt, 
"Electron  Mobility  Enhancement  by  Confining  Optical-Phonons  in  Ultrathin  GaAs/AlAs  Multiple 
Quantum  Wells,"  presented  at  the  1992  Electronic  Materials  Conference  (.June,  1992,  Cambridge, 
Mass.). 

M.  A.  Stroscio,  G.  J.  Iafrate,  K.  W.  Kim,  A.  R.  Bhatt,  M.  Dutta,  and  H.  L.  Grubin,  "Reduction  in 
Longitudinal-Optical  Phonon  Emission  Rate  in  Polar- Semiconductor  Quantum  Wires  and  Quantum 
Wells,"  presented  at  the  7th  Inti.  Conf.  on  Phonon  Scattering  in  Condensed  Matter  (August,  1992, 
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Ithaca,  New  York). 

K.  W.  Kim,  "Theory  of  Optical-Phonon  Modes  in  Ultrasmall  Electronic  Devices  (Invited),"  to  be 
presented  at  the  Advanced  Heterostructure  Transistor  Conference  (November,  1992,  Kona,  Hawaii). 
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3.  PERSONNEL 

Three  faculty  members;  Dr.  K.  W.  Kim,  Dr.  G.  J.  Iaffate,  and  Dr.  M.  A.  Stroscio,  in  the  Dept, 
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main  research  interests  include  solid-state  physics  and  electronics,  especially  the  fundamental  questions 
relating  to  the  breakdown  of  classical  solid-state  electronics  as  device  geometries  approach  the  submi¬ 
cron  and  ultrasubmicron-size  regime.  Currently,  he  is  involved  in  the  development  of  a  quantum  tran¬ 
sport  formalism  to  elucidate  tunnel-barrier  and  superlattice  hot-electron  transport  phenomena.  In  addi¬ 
tion,  he  is  also  examining  the  potential  technical  impact  of  high-temperature  superconductivity  on  elec¬ 
tronic  applications.  Dr.  M.  A.  Stroscio  has  been  associated  with  NCSU  and  with  U.S.  Army  Research 
Office  since  1985,  most  recently  as  Adjunct  Professor  of  Electrical  Engineering  and  as  Senior  Research 
Scientist,  respectively.  His  main  research  area  is  quantum  transport  theory  including  path  integral  for¬ 
malism  and  solid-state  dynamics.  He  is  currently  engaged  in  the  study  of  optical  phonon  modes  in 
heterostructures,  especially  the  effects  of  reduced  dimensionality. 

Dr.  Iafrate  and  Dr.  Stroscio  are  actively  involved  in  the  direction  of  a  post-doctoral  associate  and 
a  graduate  student  along  with  Dr.  Kim.  Due  to  the  ties  with  the  Army  Research  Office,  the  research 
activities  by  both  Dr.  Iafrate  and  Dr.  Stroscio  related  to  this  proposal  have  been  at  no  cost  to  the 
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Novel  quantum-effect  polar-semiconductor  devices  underlie  technologies 
portending  dramatic  enhancements  in  the  capability  to  process  information  orders- 
of-magnitude  faster  than  is  possible  currently.  In  many  of  these  quantum-effect 
devices,  charges  are  transported  in  quantum-well  and  quasi-one-dimensional 
quantum  wires  which  must  support  the  transport  of  charges  at  high  mobilities. 
However,  it  has  recently  been  demonstrated  that  longitudinal-optical  (LO)  phonons 
established  at  quantum-wire  interfaces  lead  to  dramatic  enhancements  in  mobility¬ 
degrading  carrier-phonon  interactions.  In  novel  mesoscopic  de  Broglie  wave 
devices,  LO  phonons  are  confined  in  structures  with  complex  geometries  which 
frequently  incorporate  metal-semiconductor  interfaces.  This  paper  outlines  the 
theory  of  confined  and  interface  LO  phonons  in  mesoscopic  devices  and 
demonstrates  that  phonon  modes  may  be  tailored  through  the  judicious  use  of 
metal-semiconductor  interlaces  in  such  a  way  as  to  dramatically  reduce  unwanted 
emission  of  interface  LO  phonons. 

Modem  techniques  for  making  semiconductor  structures  with  dimensional 
features  controlled  on  the  nanometer  scale  have  been  essential  in  leading  to 
fundamental  discoveries  as  well  as  in  stimulating  concepts  for  future  information 
processing  systems  based  on  the  exploitation  of  quantum  effects  occurring  in  such 
structures  [1-3].  As  originally  proposed  by  Sakaki  [4],  the  predicted  high  mobilities 
of  quasi-one-dimensional  wire-like  regions  of  semiconducting  material  underlie 
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Abstract 

The  dielectric  response  function  in  the  valence  bands  of  a  zincblende  semiconductor  has 
been  evaluated  in  the  random  phase  approximation  (RPA)  incorporating  the  finite  lifetime 
of  holes  at  finite  temperature.  Scattering  rates  of  minority  electrons  have  been  calculated  in 
the  Bom  approximadon.  Our  results  show  that  the  collective  morion  of  holes  (plasmons)  is 
heavily  damped  due  to  the  finite  lifetime  of  holes.  Scattering  rates  of  minority  electrons  are 
increased  due  to  these  effects. 

1.  Introduction 

An  understanding  of  the  behavior  of  electrons  in  highly  doped  semiconductors  is  increasing 
in  its  importance,  especially  for  the  design  of  high  speed  HBT’s  and  HETs,  since  electron 
transport  in  the  base  region  is  one  of  the  critical  factors  limiting  the  maximum  switching 
speed  of  such  devices.  It  is  well  known  that  the  dynamic  response  of  a  free  electron  gas  is 
described  by  the  Lindhard  dielectric  function  in  die  RPA  for  the  electron  lifetime  te+o. 
[1].  The  free  electron  gas  formula  is  easily  applicable  to  heavily  doped  n-type 
semiconductors  since  they  have  a  single  conduction  band  with  s-like  symmetry.  For  p-type 
semiconductors,  the  situation  is  much  more  complicated  because  of  the  p3/2-like  symmetry 
of  the  valence  bands.  Recently,  D.  Yevic  et  al.  evaluated  the  Lindhard  dielectric  function 
for  a  p-type  semiconductor  at  finite  temperature.  These  authors  showed  that  the  structure  of 
the  spectral  density  function  is  quite  different  from  that  for  zero  temperature.  The  scattering 
rates  of  minority  electrons  were  also  calculated  [2,  3].  However,  they  still  neglected  the 
effects  of  the  finite  hole  lifetime  even  for  finite  temperature.  Although  the  Lindhard 
dielectric  function  inherently  incorporates  the  effects  of  single  panicle  excitation  (Landau 
damping)  including  both  intra-  and  interband  transitions,  it  does  not  take  into  account  the 
finiteness  of  particle  lifetime  due  to  the  collision  effects  such  as  impurity  scattering  and 
phonon  scattering.  At  finite  temperature,  these  effects  are  very  important  since  such 
scattering  disturbs  the  collective  motion  of  quasi-panicles  and  therefore  the  dynamic 
response  function  will  be  affected  significantly.  The  purpose  of  this  paper  is  to  demonstrate 
the  effects  of  the  finite  hole  lifetime  on  the  dynamic  response  of  a  hole  gas,  as  well  as  the 
effects  on  the  lifetime  of  minority  electrons  in  p-type  semiconductors. 

2.  Theory 

From  the  fluctuation-dissipation  theorem,  the  total  inelastic  scattering  rate  for  an  electron 
with  energy  £  in  the  conduction  band  is  given  by 
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We  present  a  review  of  our  recent  results  on  electron  dynamics  in 
high  and  low  electric  fields  in  multisubband  quasi-one-dimensional 
quantum  wire  (QWI)  structures. 

First  we  derive  expressions  for  the  rates  of  electron  intra-  and 
intersubband  scattering  by  surface  optical  (SO)  and  confined 
longitudinal  optical  (LO)  phonons  in  QWIs.  Numerical  results  are 
obtained  for  a  GaAs  QWI  embedded  in  AlAs  [1,2].  There  are  two 
branches  of  SO  phonons  one  being  related  to  GaAs  and  another  to 
AlAs.  Their  frequencies  are  considerably  different  and  both  differ 
from  the  LO  phonon  frequency.  It  is  demonstrated  that  electron 
intrasubband  scattering  is  mainly  due  to  LO  phonons  whereas 
electron  intersubband  transitions  under  certain  conditions  are 
strongly  influenced  by  SO  phonons.  These  conditions  occur  when  the 
intersubband  energy  separation  in  QWIs  coincide  with  the  SO  phonon 
energy  (so  called  resonant  condition).  In  this  case  the  rate  of 
intersubband  electron  scattering  assisted  by  SO  phonons  tends  to 
infinity  for  any  electron  energy.  The  rate  of  electron-LO-phonon 
scattering  in  a  QWI  does  not  depend  crucially  on  the  intersubband 
energy  separation.  Note  that  only  "forward"  electron  scattering  can 
be  in  the  resonance.  The  "backward"  scattering  rate  depends  weakly 
upon  resonant  conditions.  High  intersubband  scattering  rates  under 
resonant  conditions  must  yield  fast  electron  redistribution  between 
subbands  in  resonance.  However,  there  is  no  momentum  relaxation 
due  to  "forward"  resonant  intersubband  electron  scattering. 
Therefore,  resonant  intersubband  scattering  should  not  affect 
electron  transport  unless  the  subbands  in  resonance 
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As  device  dimensions  in  nanoscale  structures  and  mesoscopic  systems 
are  reduced,  the  characteristics  and  interactions  of  dimensionally-confined 
longitudinal-optical  (LO)  phonons  deviate  substantially  from  those  of  bulk  polar 
semiconductors.  This  account  emphasizes  LO-phonon  effects  arising  in  three 
separate  systems:  Short-period  AJAs-GaAs  superlattices;  rectangular  GaAs 
quantum  wires  embedded  in  AlAs;  and  metal-semiconductor  interfaces  such  as 
those  in  de  Broglie  wave  quantum-interference  devices. 

I.  INTRODUCTION 

In  nanoscale  and  mesoscopic  systems,  the  effects  of  confinement  on 
carriers  have  been  studied  extensively.  However,  to  properly  model  carrier 
energy  loss  in  nanoscale  and  mesoscopic  systems,  it  is  essential  that  calculations 
of  carrier  scattering  by  longitudinal-optical  phonons  take  into  account  the  fact 
that  confinement  also  changes  the  strength  and  spatial  properties  of  longitudinal- 
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Abstract,  Macroscopic  dielectric  continuum  models  of  optical-phonon  modes 
predict  enhancements  in  the  magnitudes  of  the  surface-optical  (so)  modes  in 
double-barrier  heterostructures  as  the  heterojunction-to-heterojunction  separation 
is  reduced.  In  this  paper,  the  ratio  of  electron  scattering  by  the  so-phonon  modes 
to  that  by  the  (electrostatic)  confined  longitudinal -optical-  (LO-)  phonon  modes  is 
calculated  for  a  GaAs/AIAs  short-period  superlattice  based  on  the  assumption 
that  electron-so-phonon  scattering  may  be  described  by  a  scalar  potential.  The 
scaling  of  the  ratio  of  electron-so-phonon  scattering  to  electron-LO-phonon 
scattering  as  a  function  of  the  superlattice  period  provides  a  sensitive  test  of  the 
appropriateness  of  the  scalar-potential  model.  The  effect  of  phonon  confinement 
on  electron-optical-phonon  scattering  rates  is  presented  for  rectangular  quantum 
wires  as  well.  A  major  conclusion  of  these  new  results  is  that  it  is  essential  to 
model  phonon  confinement  properly  in  predicting  carrier  transport  properties  in 
mesoscale  structures. 


Interactions  between  carriers  and  longitudinal-optical 
(lo)  phonon  modes  in  heterostructures  are  strongly 
affected  by  the  changes  in  the  Frohlich  Hamiltonian 
caused  by  phonon  confinement  and  localization,  as  well 
as  by  the  changes  in  the  electronic  wave  function  due  to 
the  confining  potential.  The  presence  of  heterointerfaces 
gives  rise  to  the  confinement  of  optical  phonons  in  each 
layer  (i.e.,  confined  modes)  and  the  localization  in  the 
vicinity  of  interfaces  (i.e.,  interface  modes  or  surface- 
optical  (SO)  modes).  There  have  been  suggestions  that 
interactions  by  the  interface  modes  can  be  significant  and 
that  the  scattering  rate  due  to  the  confined  modes  can  be 
considerably  reduced  in  some  structures  compared  to  the 
bulk  LO-phonon  scattering  rate  [1-3].  Therefore,  an 
appropriate  treatment  of  the  optical-phonon  modes  in 
quantized  systems  is  essential  for  the  understanding  of 
electron  transport  in  heterostructures.  Recently,  both 
macroscopic  and  microscopic  approaches  to  electron- 
optical-phonon  interactions  in  heterostructures  [1-11] 
have  been  applied  in  theoretical  treatments.  Enhanced 
electron-so-phonon  scattering  in  polar  semiconductors 
with  confining  dimensions  less  than  about  50  A  has  been 
indicated  recently  [12, 13]  on  the  basis  of  scalar-potential 
modes.  The  appropriateness  of  using  a  scalar  potential  to 
model  the  electron-so-phonon  interaction  has  been 
questioned  [14].  In  this  paper,  we  model  the  electron-so- 
phonon  interaction  with  a  scalar  potential  to  calculate 
the  ratio  of  electron-so-phonon  scattering  to  electron- 
confined-LO-phonon  scattering  in  a  short-period 


GaAs/AIAs  superlattice;  the  results  scale  inversely  with 
the  superlattice  quantum-well  dimension  in  qualitative 
agreement  with  recent  experiments  [15].  Such  com¬ 
parisons  provide  convincing  evidence  that  the  so  modes 
present  in  these  experiments  contain  an  active  lo  compo¬ 
nent.  This  scalar-potential  description  is  also  applied  to 
calculate  the  electron-optical-phonon  scattering  rates  in 
GaAs  rectangular  quantum  wires. 

The  GaAs/AIAs  short-period  superlattice  considered 
in  this  study  has  the  GaAs  layer  width  a  and  the  AlAs 
layer  width  b.  The  ratio  of  electron-so-phonon  scatter¬ 
ing  to  electron-LO-phonon  scattering  in  this  superlattice 
is  calculated  through  the  Fermi  golden  rule  by  treating 
the  electron-LO-phonon  and  electron-so-phonon  inter¬ 
action  Hamiltonians  [4,6]  as  perturbation  Hamiltonians 
and  by  taking  the  electronic  wavefunction  along  the 
growth  direction  (i.e.,  the  z-direction)  as  that  given  by 
Cho  and  Prucnal  for  the  maximum  edge  of  the  first 
subband  [16].  To  obtain  an  order-of-magnitude  estimate 
of  the  ratio  of  electron-so-phonon  scattering  to 
electron-LO-phonon  scattering  in  the  GaAs/AIAs  super¬ 
lattice,  we  consider  the  ratio  of  principal  factors,  Y,  in  the 
integrands  of  the  respective  scattering  rates.  In  this 
calculation,  the  contribution  by  the  S  —  mode  is  neglec¬ 
ted  due  to  its  relatively  small  contribution  as  compared 
with  the  S+  mode  [6].  Figure  1  shows  the  ratio,  T 
plotted  as  a  function  of  X  (  =  qa  where  q  is  the  phonon 
momentum  parallel  to  the  heterointerface)  for  a 
GaAs/AIAs  short-period  superlattice  with  b  =  a  and 
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Abstract.  The  Monte  Carlo  simulation  of  electron  transport  in  multisubband 
quasi-one-dimensional  GaAs/AIAs  quantum  wires  (qwis)  are  presented.  The 
electron  intrasubband  and  intersubband  scattering  by  surface-optical  and 
confined  longitudinal-optical  phonons  in  qwis  has  been  included  in  the  program. 
It  is  demonstrated  that  at  room  temperature  the  electron  drift  velocity  in  the  QWi  is 
suppressed  by  electron  intersubband  scattering  and  does  not  exceed  the  bulk 
material  values.  The  energy  dependence  of  the  total  scattering  rate  in  ideal  qwis 
exhibits  multiple  sharp  peaks  related  to  intersubband  transitions.  The  scattering 
rates  in  the  real  qwis  with  variable  thicknesses  are  calculated  as  well.  The  results 
show  that  even  small  variation  in  thickness  leads  to  a  significant  broadening  of 
the  very  first  peaks  and  complete  washing-out  of  the  peak-like  structure  at  higher 
energies. 


1.  Introduction 

During  the  years  following  Sakaki’s  work  [1],  a  number 
of  papers  on  quantum  wires  (qwis)  appeared,  dealing 
mainly  with  theoretical  analyses  of  quantization,  scatter¬ 
ing  mechanisms  and  linear  transport  in  low  electric  fields 
(see  e.g.,  [2]  and  papers  in  [3]).  Most  of  the  previous 
studies  on  nonlinear  electron  transport  in  qwis  have 
been  limited  to  the  case  of  scattering  by  bulk  three- 
dimensional  (3D)  phonon  modes  [4],  However,  experi¬ 
ments  clearly  demonstrate  the  importance  of  surface- 
optical  (so)  modes  [5]  and  phonon  confinement  [6]  in 
qwis.  Moreover,  most  theoretical  studies  on  qwis  deal 
with  ideal  ID  systems  characterized  by  fixed  subband 
energy  positions.  However,  all  current  and  foreseeable 
future  technologies  of  fabricating  qwis  do  not  assure  the 
possibility  to  create  ideal  structures  with  constant  thick¬ 
ness  [7-10].  The  variation  of  qwi  thickness  results  in  the 
variation  of  subband  energy  positions.  Consequently, 
electron  scattering  is  no  longer  energetically  coherent  in 
different  parts  of  a  qwi  and  this  should  lead  to  the 
broadening  or  even  complete  washing-out  of  the  reso¬ 
nant  peaks  [8-10], 

In  the  present  paper  we  consider  nonlinear  (hot) 
electron  transport  in  a  rectangular  GaAs  qwi  embedded 
in  AlAs.  Both  the  phonon  confinement  and  multi¬ 
subband  structure  are  taken  into  account.  We  consider 
electron  scattering  in  an  ideal  qwi  with  a  constant 


thickness  along  the  wire  as  well  as  in  a  real  qwi  with 
variable  thickness.  The  Monte  Carlo  simulations  are 
performed  for  the  ideal  qwi. 


2.  Scattering  rates 

The  model  of  the  qwi  and  the  scattering  rates  is  discussed 
in  detail  in  [10. 1 1  ].  We  consider  ideal  qwis  with  constant 
thickness  and  real  qwis  with  a  smooth  variation  of 
thickness  along  the  structure  so  that  the  characteristic 
length  of  the  fluctuations  is  much  greater  than  the  de 
Broglie  wavelength  AB.  Numerical  calculations  of  the 
scattering  rates  have  been  performed  for  the  GaAs  qwi  of 
dimensions  Ly  =  150  and  L.  =  250  A  embedded  in  AlAs. 
We  have  considered  9  subbands.  The  scattering  rates  in 
real  qwis  have  been  calculated  assuming  a  thickness 
varying  with  the  harmonic  law.  The  variation  amplitudes 
5Ly  and  5L.  have  been  chosen  from  0.01  to  0.1  fraction  of 
the  corresponding  thicknesses  Ly  and  L..  We  have 
considered  three  cases  with  different  mutual  phases  of  the 
variation  of  Ly  and  L.:  L.  =  L.0  +  SL.coMKx), 
L.  =  L:0  +  &L.sm{Kx)  and  L.  =  L:0  +  5L.cos(2Kx) 
with  Ly  =  Lyo  +  5Lycos(Kx)  in  all  cases.  It  has  been 
found  that  the  results  do  not  show  significant  difference 
between  these  models  as  long  as  we  are  in  the  iiir.it  of 
smooth  fluctuations  KAB«1.  The  calculated  energy 
dependences  of  the  total  electron  transition  rate  from  the 
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Abstract.  The  role  of  confined  phonon  modes  in  determining  the  energy 
relaxation  of  hot  electrons  in  low-dimensional  semiconductor  microstructures  is 
discussed  within  a  dielectric  continuum  model  for  the  lo  phonon  confinement 
and  a  long  wavelength  Frohlich  model  for  the  electron -phonon  interaction. 
Numerical  results  are  provided  for  the  hot-electron  relaxation  rate  as  a  function  of 
electron  temperature  and  density  for  GaAs  quantum  wells  and  quantum  wires  by 
taking  into  account  emission  of  slab  phonon  modes.  Comparison  with  existing 
experimental  results  shows  some  evidence  for  slab  phonon  emission  in  inter¬ 
subband  electronic  relaxation  in  reasonably  narrow  quantum  wells.  It  is  argued 
that  most  experiments  can  be  interpreted  in  terms  of  an  electron-bulk  phonon 
interaction  model  (i.e.  by  taking  into  account  the  effect  of  confinement  only  on 
the  electrons  and  assuming  the  phonons  to  be  the  usual  bulk  three-dimensional 
phonons)  because  a  number  of  important  physical  processes,  such  as  screening, 
the  hot  phonon  effect,  phonon  self-energy  correction  etc,  make  it  difficult  to 
distinguish  quantitatively  between  various  models  for  phonon  confinement, 
except  perhaps  in  the  narrowest  (<50A)  wells  and  wires.  Detailed  numerical 
results  for  the  calculated  intra-subband  relaxation  rate  in  GaAs  quantum  wires  are 
provided  within  the  slab  phonon  and  the  electron  temperature  model,  including 
the  effects  of  dynamical  screening,  quantum  degeneracy  and  non-equilibrium  hot 
phonons. 


1.  Introduction 

A  large  number  of  theoretical  papers  have  recently 
appeared  [1-21]  on  the  subject  of  interaction  between 
free  electrons  and  optical  phonons  in  confined  low¬ 
dimensional  semiconductor  structures  such  as  semicon¬ 
ductor  heterostructures,  quantum  wells,  superlattices 
and  quantum  wires.  While  it  is  well  established  [1]  that 
the  effect  of  quantum  confinement  on  the  electronic 
degrees  of  freedom  must  be  included  in  any  reasonable 
theoretical  model  of  electron-phonon  interactions  in 
these  microstructures,  the  recent  emphasis  in  the  litera¬ 
ture  has  been  on  the  inclusion  of  confinement  effects  on 
the  phonon  modes  of  the  structures  and  the  consequent 
modification  of  the  Frolich  Hamiltonian  describing  the 
electron-phonon  interaction  in  the  system.  Because 
Frohlich  interaction  plays  an  important  quantitative  role 
in  determining  many  aspects  of  the  electronic  properties 
(e.g.  electron  mobility,  hot-electron  relaxation,  inelastic 
mean  free  path,  polaron  effective  mass)  of  micro- 
structures,  the  possible  influence  of  confined  phonon 
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modes  in  semiconductor  microstructures  is  of  some 
interest.  The  purpose  of  this  paper  is  to  discuss  the 
quantitative  role  of  confined  phonon  modes  from  the 
perspective  of  one  particular  class  of  phenomena, 
namely,  the  hot-electron  energy  relaxation  process  in 
quantum  wells  and  wires.  Our  main  concern  in  this 
article  is  the  role  of  confined  phonon  modes  in  the  energy 
relaxation  of  hot  electrons  in  GaAs  microstructures  with 
the  particular  emphasis  on  qualitative  and  quantitative 
understanding  of  recent  experimental  results  [22-30]. 

The  most  unambiguous  evidence  for  the  existence  of 
confined  phonon  modes  in  microstructures  comes  from 
Raman  scattering  experiments  [22-24]  in  thin  GaAs- 
AlAs  superlattices  where  confinement  effects  on  both  the 
acoustic  phonons  [22]  and  optical  phonons  [23, 24]  have 
directly  been  observed  for  layer  thicknesses  typically 
below  50  A.  Unfortunately,  such  Raman  measurements 
cannot  tell  us  much  about  the  role  of  confined  phonon 
modes  in  determining  various  electronic  properties.  In 
addition,  the  Raman  measurements  (and  this  limitation 
applies  also  to  all  of  the  ‘so-called’  first  principles  or 
microscopic  lattice-dynamical  calculations)  are  neces¬ 
sarily  on  superlattices  whereas  the  electronic  experiments 


Oji  ■o'*  ' 


Philosophical  Magazine  Letters,  1992,  Vol.  65,  No.  4,  173-176 


Dramatic  reduction  in  the  longitudinal-optical  phonon  emission  rate 
in  polar-semiconductor  quantum  wires 

By  Michael  A.  StroscioI,  K.  W.  KimJ,  Gerald  J.  Uprated 
Mitra  Dutta§  and  Harold  L.  Grubin|| 

fU.S.  Army  Research  Office,  P.O.  Box  12211, 

Research  Triangle  Park,  North  Carolina  27709,  U.S.A. 
t  Department  of  Electrical  and  Computer  Engineering, 

North  Carolina  State  University,  Raleigh,  North  Carolina  27695,  U.S.A. 

§U.S.  Army  Electronics  Technology  and  Devices  Laboratory, 

Ft  Monmouth,  New  Jersey  07703,  U.S.A. 

||  Scientific  Research  Associates,  Inc.,  50  Nye  Road, 

Glastonbury,  Connecticut  06033,  U.S.A. 

[Received  8  November  1991  and  accepted  17  December  1991 ] 

Abstract 

Novel  quantum-effect  polar-semiconductor  structures  underlie  technologies 
portending  dramatic  enhancements  in  the  capability  to  process  information  orders 
of  magnitude  faster  than  is  possible  currently.  In  many  embodiments  of  these 
quantum-effect  structures,  charges  are  transported  in  quasi-one-dimensional 
quantum  wires  which  must  support  the  transport  of  charges  at  high  mobilities. 
However,  it  has  recently  been  demonstrated  that  the  longitudinal-optical  (LO) 
phonons  established  at  quantum-wire  interfaces  lead  to  dramatic  enhancements  in 
carrier-phonon  interactions  and  concomitant  degradation  in  carrier  mobility.  This 
letter  demonstrates  that  phonon  modes  may  be  tailored  through  the  judicious  use  of 
metal-semiconductor  interfaces  in  such  a  way  as  to  dramatically  reduce  unwanted 
emission  of  interface  LO  phonons  and,  consequently,  to  lead  to  the  achievement  of 
high  quantum-wire  mobility. 


Modern  fabrication  techniques  for  making  semiconductor  structures  with 
nanometre-scale  dimensional  features  have  been  essential  in  leading  to  fundamental 
discoveries  (Esaki  1974,  Stein,  von  Klitzing  and  Weimann  1983)  as  well  as  in 
stimulating  concepts  for  future  information-processing  systems  based  on  the  exploit¬ 
ation  of  quantum  effects  occurring  in  such  structures  (Capasso  and  Datta  1990).  As 
originally  proposed  by  Sakaki  (1980),  the  predicted  high  mobilities  of  quasi-one- 
dimensional  wire-like  regions  of  a  semiconducting  material  underlie  many  proposed 
quantum-wire  system  concepts  (Luryi  and  Capasso  1985,  Sakaki  1989)  and  have 
resulted  in  pioneering  efforts  leading  to  the  actual  fabrication  of  quantum  wires  (Watt, 
Sotomayor  Torres,  Arnot  and  Beaumont  1990).  Furthermore,  quantum-wire  struc¬ 
tures  have  been  essential  to  seminal  studies  underlying  quantum-wire  laser  concepts 
(Tsuchiya  et  al.  1989)  and  quantum-coupled  electronic  systems  (Reed  et  al.  1988). 
Recently,  however,  theoretical  studies  of  the  interaction  between  longitudinal-optical 
(LO)  phonons  and  carriers  in  a  polar-semiconductor  quantum  wire  (Stroscio  1989) 
have  revealed  the  presence  of  discrete  LO  phonon  modes  similar  to  those  identified  for 
polar-semiconductor  quantum  wells  (Kliewer  and  Fuchs  1966,  Licari  and  Evrard  1977, 
Mori  and  Ando  1989,  Kim  and  Stroscio  1990).  As  for  the  case  of  quantum  wells, 
interface  LO  phonons  are  established  at  the  semiconductor-semiconductor  bound¬ 
aries  of  quantum  wires  (Kim  et  al.  1991).  In  addition,  it  has  been  shown  that,  for  carrier 
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The  rates  of  electron  intra-  and  intersubband  scattering  by  surface  optical  (SO)  and 
confined  longitudinal  optical  (LO)  phonons  in  quasi-one-dimensional  GaAs/AlAs 
q  tin  turn  wires  (QWIs)  are  calculated.  It  is  shown  that  electron-SO-phonon  intersubband 
scattering  can  be  resonant,  so  that  the  scattering  rate  tends  to  infinity  when  intersubband 
energy  separation  approaches  SO  phonon  energy.  The  energy  dependence  of  the  total 
scattering  rate  in  ideal  QWIs  exhibits  multiple  sharp  peaks  related  to  intersubband 
transitions.  The  scattering  rates  in  the  real  QWIs  with  variable  thickness  are  calculated. 
The  results  show  that  even  small  variation  in  thickness  leads  to  the  significant  broadening 
of  the  very  first  peaks  and  complete  washing-out  of  the  peak-like  structure  at  higher 
energies.  Monte  Carlo  simulations  of  electron  transport  in  the  QWI  have  been  performed. 
It  is  demonstrated  that  electron  drift  velocity  in  the  QWI  is  considerably  suppressed  by 
electron  intersubband  scattering  and  is  considerably  lower  than  the  bulk  material  values. 
The  role  of  SO  phonons  in  electron  energy  dissipation  is  discussed. 


I.  Introduction 

Despite  rapidly  growing  number  of  publications  on 
quantum  wires  (QWIs),  the  theoretical  investigations  are 
limited  either  to  the  case  of  extreme  quantum  limit  (EQL) 
wherein  only  one  subband  is  considered  (1,21  or  (and)  to 
the  case  of  scattering  by  bulk  three-dimensional  (3-D) 
phonon  modes  [3].  However,  it  is  difficult  to  meet  these 
limitations  because  elections  can  populate  upper  subbands 
at  higher  temperatures  or  in  the  hot-electron  regime  (4)  and 
experiments  evidently  demonstrate  importance  of 
surface-optical  (SO)  modes  [5]  and  phonon  confinement 
[6J  in  QWIs.  Moreover,  most  theoretical  studies  on  QWIs 
deal  with  ideal  1-D  systems  characteriz'd  by  fixed  subband 
energetic  positions.  The  unique  feature  of  ideal  QW 
the  well-pronounced  resonant  nature  of  electron  scattc  n.. 
as  a  result  of  multiple  sharp  peaks  (diverging  to  infinity) 
on  energy  dependence  of  the  total  electron  scattering  rate. 
However,  all  current  and  foreseeable  future  technologies 
of  fabricating  QWIs  do  not  assure  the  possibility  to  create 
ideal  structures  with  constant  thickness  [4,7].  The 
variation  of  QWI  thickness  results  in  the  variation  of 
subband  energetic  position.  Consequently,  electron 
scattering  is  no  longer  energetically  coherent  in  different 
pans  of  a  QWI  and  this  should  lead  to  the  broadening  or 
even  complete  washing-out  of  the  resonant  peaks  (8,9]. 

The  aim  of  the  present  paper  is  to  investigate  both  the 
phonon  confinement  and  multi-subband  structure  in  a 
rectangular  QWI  of  polar  semiconductor  embedded  in 
another  polar  semiconductor.  We  consider  an  ideal  QWI 


with  thickness  constant  along  the  wire  (Section  II)  and  a 
real  QW  with  variable  thickness  (Section  III).  In  Section 
IV  some  preliminary  results  of  Monte  Carlo  simulation  are 
discussed 


II.  Ideal  QWI 

The  1  -D  electron  energy  in  a  rectangular  QWI  is  of 
well-known  form 


E(q,j,l)  =  E,(qx)  +  Eji , 


where  E,  is  the  electron  kinetic  energy,  q,  being  the 
electron  wave-number  in  x-direction,  and  Ejt  is  the 
subband  (j J }  energetic  position  with  respect  to  the  bulk 
ground  level.  It  is  seen  from  Eq.(l)  that  subband  energetic 
positions  depend  on  the  structure  thickness  L^  and  L,.  In 
this  Section  we  will  consider  an  ideal  QWI  with  constant 
thickness  along  the  wire,  so  that  E^,  is  independent  of  x. 
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Abstract.  The  present  study  deals  with  electron  intersubband  scattering  in  teal  quantum 
wire  structures.  Both  the  mulli-subband  structure  and  confined  phonon  modes  arc 
considered  together.  The  rates  of  scattering  by  confined  longitudinal-optical  (tx>)  phonons 
and  by  surface-optical  (SO)  phonons  are  calculated  taking  into  account  all  possible  lo 
phonon  modes  as  well  as  all  possible  electron  intersubband  transitions.  The  estimations  of 
transition  rates  for  GaAs/AlAs  owls  have  shown  that  inlrasubband  electron  scattering  and 
most  intersubband  transitions  are  due  primarily  to  scattering  by  confined  io  phonons,  but 
in  resonant  intersubband  transitions  the  contribution  of  so  phonons  may  be  dominant 
when  the  phonon  energy  is  close  to  the  intersubband  energy  separation.  Moreover, 
dectron-so-phonon  scattering  might  play  an  important  pan  in  low-temperature  dectron 
transport  because  the  GuAs-like  so  mode  is  shifted  towards  lower  frequencies  compared 
with  that  of  lo  phonons.  The  energy  dependence  of  the  total  scattering  rate  in  an 
ideal  quantum  wire  exhibits  multiple  sharp  peaks  related  to  each  intersubband  transition. 
These  peaks  originate  from  the  resonant  nature  of  the  density  of  states  in  ideal  one¬ 
dimensional  rystems.  It  is  demonstrated  that  in  real  quantum  wires  with  variable  thickness 
the  resonant  peaks  broaden  or  even  disappear  due  to  variation  of  tubband  energies. 


I.  Introduction 

The  progress  in  semiconductor  technology  has  provided  the  means  to  fabricate  the  so- 
called  quantum  wires  (qwis)  with  quasi-one-dimensional  (id)  structures.  It  has  been 
suggested  that  qwis  will  exhibit  carrier  mobilities  well  above  IQ6  cm2  V-1  s"1  [1], 
but  these  high  values  of  the  mobility  have  not  yet  been  observed  experimentally.  The 
expected  enhancement  of  the  carrier  mobility  in  qwis  should  stem  from  the  restriction 
of  momentum  space  to  one  dimension  as  well  as  the  resulting  reduction  of  final  states 
for  scattered  electrons.  This  point,  however,  needs  to  be  clarified.  Despite  the  rapidly 
growing  number  of  publications  on  qwis,  the  theoretical  investigations  of  electron 
transport  controlled  by  optical  phonon  scattering  are  limited  either  to  the  case  of  the 
extreme  quantum  limit  (EQL)  wherein  only  one  subband  is  considered  (2-4}  and/or  to 
the  case  of  scattering  by  bulk  three-dimensional  (3D)  phonon  modes  [5—7].  However, 
due  to  technological  limitations  the  confinement  of  electrons  is  relatively  weak  and 
electrons  can  populate  upper  subbands  at  higher  temperatures  or  in  the  hot-electron 
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In  this  paper,  the  commonly  used  but  idealistic  formulation  of  quantized  optical-phonon  modes 
for  a  heterostructure  system  with  only  two  heterojunctions  (i.e.,  single  quantum-well  structures) 
is  extended  to  the  more  realistic  case  of  multiheterointerface  structures.  By  applying  the 
macroscopic  dielectric  continuum  approach,  dispersion  relations  and  interaction  Hamiltonians 
for  interface-phonon  modes  are  derived  for  a  double-barrier  structure  and  scattering  rates  based 
on  these  results  are  used  to  determine  the  range  of  practical  validity  of  the  idealistic  model  using 
interaction  Hamiltonians  appropriate  for  single  quantum  wells  with  infinite  barrier  widths.  It  is 
found  that  when  the  dimensions  of  the  structures  are  larger  than  approximately  30  A,  this 
simplified  description  can  be  applied  to  multiheterointerface  structures  in  general  with 
reasonable  accuracy. 


I.  INTRODUCTION 

It  is  well  known  that  electron  interactions  with 
longitudinal-optical  (LO)  phonon  modes  in  heterostruc¬ 
tures  are  strongly  affected  by  the  changes  in  the  Frohlich 
Hamiltonian  caused  by  phonon  confinement  and  localiza¬ 
tion,  as  well  as  by  the  changes  in  the  electronic  wave  func¬ 
tion  due  to  the  confining  potential.  The  presence  of  hetero- 
interfaccs  gives  rise  to  the  confinement  of  optical  phonons 
in  each  layer  (i.e.,  confined  modes)  and  the  localization  in 
the  vicinity  of  interfaces  (i.e.,  interface  modes).  There  have 
been  suggestions  that  carrier  interactions  with  the  interface 
modes  can  be  significant  and  that  the  scattering  rate  due  to 
the  confined  modes  can  be  considerably  reduced  in  some 
structures  compared  to  the  bulk  LO-phonon  scattering 
rate.1-5  Therefore,  an  appropriate  treatment  of  the  optical- 
phonon  modes  in  quantized  systems  is  essential  for  the 
understanding  of  electron  transport  in  heterostructures. 
Recently,  both  macroscopic1'10  and  microscopic1112  ap¬ 
proaches  of  electron-optical-phonon  interactions  in  hetero- 
structures  have  been  applied  in  theoretical  treatments.  De¬ 
tailed  microscopic  calculations  of  optical-phonon  modes  in 
polar  semiconductors  indicate  that  the  dielectric  contin¬ 
uum  model  provides  an  accurate  formalism  for  modeling 
electron-optical-phonon  interactions.12  However,  most  of 
these  theoretical  analyses  have  been  confined  to  highly 
symmetric  and/or  simple  structures  such  as  single  quan¬ 
tum  wells  or  superlattices.  Application  of  even  a  simple 
macroscopic  model,  not  to  mention  the  microscopic  ab 
initio  models,  becomes  highly  complicated  due  to  the  cou¬ 
pling  between  adjacent  interfaces  when  the  structure  has 
multiple  heterointerfaces  or  is  asymmetric.  There  are  many 
such  structures  in  which  a  rather  accurate  description  of 
the  electron-optical-phonon  interaction  (such  as  the  role  of 
interface  phonons)  is  of  prime  concern.11'18  Thus,  it  is 
desirable  to  investigate  the  possibility  of  utilizing  an  ap¬ 


proximate  description  of  phonon  modes  and  related  inter¬ 
action  Hamiltonians,  while  retaining  the  quantized  nature 
of  phonon  modes  in  the  multiheterointerface  structures. 

In  this  paper,  we  report  a  detailed  study  of  the 
electron-optical-phonon  interaction  in  a  double-barrier 
heterostructure  to  explore  the  possibility  of  achieving  a 
simplified  description  such  as  discussed  above.  The  double¬ 
barrier  structure  has  been  chosen  as  a  specific  example  due 
to  its  structural  similarity  (with  an  added  complexity  of 
finite  barrier  width)  to  single  quantum  wells  (i.e.,  infinite 
barrier  width)  for  which  quantized  optical-phonon  modes 
are  already  well  known,  and  its  prominence  in  recent  stud¬ 
ies  which  demonstrated  the  existence  and  importance  of 
interface-phonon  modes.11  16  In  this  study,  the  interaction 
Hamiltonians  are  derived  by  applying  the  macroscopic  di¬ 
electric  continuum  model4-10  without  any  further  approx¬ 
imations.  The  resulting  scattering  rates  are  compared  with 
those  obtained  by  assuming  the  Hamiltonians  appropriate 
for  single  quantum  wells  for  which  finite  barrier  widths  are 
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FIG  I  Schematic  drawing  of  the  double-barrier  heterostructure  dis¬ 
cussed  in  this  study  The  z  axis  is  chosen  as  the  growth  direction  of  the 
lattice. 
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